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Abstract. Polarization measures, on wide angular scales, together with anisotropy data, can 
fix DE parameters. Here we discuss the sensitivity needed to provide significant limits. Our 
analysis puts in evidence that a class of models predicts low correlation or anticorrelation between 
'• polarization and anisotropy at low I. This class includes open models and models with DE due 

to a Ratra-Peebles (RP) potential. Results on this point, given in a previous paper of ours, 
are updated and partially corrected. We outline that, with the sensitivity of experiments like 
SPOrt or WMAP, high values of A (energy scale in the RP potential) can be excluded. With 
Q-r the sensitivity expected for PLANCK, the selection will extend to much lower A's. 
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1 Introduction 

r—i ■ 

The nature of Dark Energy (DE) is one of the main puzzles of cosmology. DE was first required by 
£T) | SNIa data [1], but a flat Universe with Q m ~ 0.3 and O^/i 2 ~ 0.02 is also favored by CBR and LSS 
observations [2,3] (f2 m &: matter, baryon density parameters; h: Hubble parameter in units of 100 
km/s/Mpc; CBR: cosmic background radiation; LSS: large scale structure). 

DE could be a false vacuum; then, its pressure and energy density (pde and pde) have ratio 
w = — 1. This however requires a severe fine tuning at the end of the EW transition. Otherwise, DE 
can be a scalar field <fi self-interacting through a potential V(4>) (dynamical DE [4,5,6]). Then 

P< pde = <j) 2 /2a 2 + V{4>) , p DE = 4?/2a? -V{(f) (1.1) 

(derivatives are in respect to the conformal time t). As soon as pk = 4> 2 /2a 2 < V, it is w < 
For pk/V ~ 1/2, it is w ~ —1/3 and dynamical DE approaches an open CDM behavior. Smaller 
Pk/V ratios approach w = —1 and a ACDM behavior. To work out w(a), the Friedman equations, 
together with the equation of <j), are to be integrated; the solutions depend on the shape of V, which, 
■ in principle, is largely arbitrary. 

Among potentials admitting a tracker solution, the RP [5] and SUGRA [6] expressions 

. 5? ! V (cf>) = A 4 + a /4> a , V{4>) = (A 4+c 7 4> a ) exp(4^ 2 /™') , (1.2) 

are particularly relevant, as they originate within the frame of Supersymmetric (SUSY) theories. Here, 
A is an energy scale, set in the interval 1-10 12 GeV; m p is the Planck mass. Once A and £Ide are 
fixed, the exponent a is set. RP and SUGRA potentials yield fast and slowly varying w, respectively. 

Dynamical DE and ACDM often predict similar observational outputs. This is welcome, as ACDM 
is a good fit to data. However, here we show that measures of anisotropy and polarization of CBR, 
at large angular scales, constrain V and distinguish dynamical DE from ACDM. Accordingly, exper- 
iments in progress already exclude some parameter range for RP models and higher sensitivities can 
discriminate even better. In a previous paper [7], similar results were provided; that paper, however, 
contained a numerical mistake. Here we correct some of its quantitative results. 

Our procedure includes a likelihood analysis, assuming polarization data provided by the Sky 
Polarization Observatory [8], both with the expected experimental noise level and a higher sensitivity. 
Results, however, can be straightforwardly extended to other observational contexts. 



2 CMB angular spectra from the Boltzmann equations; theory 

T E TE 

The angular CBR spectra C l ' ' (only E-mode is considered through this paper) can be worked 
out from the linear fluctuation evolution, obtained by solving the Boltzmann equation for the photon 



distribution. The treatment is discussed in a number of papers [9], giving also the equations for the 
other components of a model. In this section we show that these equations, in several cases, yield a 
low or negative Cf E for low I. All definitions used are the same as in CMBFAST [10]. 

This effect arises because of the simultaneous action of the ISW effect and of the opacity r = 
$1° n e (t')a T a(t') dt' '. (Notice that — f = an e aT)', n- e and <jt are the free electron density and the 
Thomson cross-section. ISW effect arises when we pass either from matter to curvature dominance 
(open models) or from matter to vacuum dominance (ACDM models). However, only in the former 
case and in the presence of opacity, anticorrelation arises. Some RP models induce anticorrelation 
because their features more closely approach open CDM, rather than ACDM. 

Let us then indicate by Fi(k,t) and Gi(k,t) the Boltzmann components for anisotropy and polar- 
ization (k is the wave-number). The equations for the Gi components in flat models read: 

G l = -f[ — Gi + H (5|0 + S f)] + _A_[/ Gi-i -{1 + l)G l+1 ] (2.1/) 

(8i n is the Kronecker symbol). Here 

U = G Q + G 2 + F 2 , 

is the only vehicle from anisotropy to polarization. In open models, eq. (2. If) becomes: 

Gi = -t[ -Gi + 5(^ + 5 -f)} + ^^[Ibi G,_! -(1 + 1)6 J+1 G, +1 ] ; (2.1o) 

here (3 2 = k 2 + K and bf = 1 — Kl 2 /(3 2 , with K = -(1 — Q m )HQ, H is today's Hubble parameter. 

Initially all Gi terms are zero; to switch them on, the quadrupole F 2 (k,t) must be great when n e 
is not so low. For wavelengths 2ir/k entering the horizon well after recombination, F 2 (k,t) switch on 
when n e has almost vanished, unless reionization occurs. Notice that the horizon size at recombination 
corresponds to I <C 200. Without reionization, below such I, we expect low Gi's. 

Let F° (Gf ) be the present value of harmonics. In spatially flat models, the angular spectra read 

C f = \j d 3 kP Q (k) \F?(k)\ 2 ,Cf = ^j d 3 kP Q (k) \GUk)\ 2 , Cf E = ^J d 3 kP (k)Fnk)GUk), 

(2.3/) 

P (k) being the primeval fluctuation spectrum. In open model, instead, they read 

Cf = \J d 3 (3P (q) \F?m 2 ,Cf = ^J d 3 (3P (q) \Gf((3)\ 2 , Cf E = ^f d 3 (3P (q) F° ((3)G° ((3) ; 

(2.3o) 

here q = (f3 2 - 4K) 2 //3(/3 2 - K). Clearly, for Q m = 1, both (3 and q return k. 

Comparing eqs. (2.1) and (2.3) shows suitable shifts in the A:-space. In particular, the bi coefficients 
cause a shift of C/ peaks, while the passage from k to (3 and q, in eq. (2.3o), displace the power through 
the harmonics Fi at small /, i.e. on scales comparable with the curvature scale. Apart of these shifts, 
the gravitational field fluctuations (h and rj) obey similar but different equations: 

Open models : 2k 2 fj = 87rGa 2 \^2(p c + p c )9 c — hp 0iCr (l — Q m )/a 2 ]. (2-4o) 

c 

Flat models : 2k 2 i] = ^Ga 2 ^{ Pc + Pc )9 c + (pde + Pde)0de}- (2.4/) 

c 

Here Yj C sums over an (relativistic or non-relativistic) matter components apart of DE; p 0)Cr is today's 
critical density. The wave-numbers k 2 and k 2 are just shifted by 3K. 

Notice that, in a ACDM model, no DE fluctuations exist while pde = ~Pde, so that the second 
term at the r.h.s. of eq. (2.4f) vanishes. In models with dynamical DE, instead, Ode 7^ and 
Pde 7^ —Pde- Accordingly, the second term in square brackets in eq. (2.4f) may read 6dePo,ct{^ — 
^m)(l + w){pDE / Po,de) (the last parenthesis tells us how DE energy scales with a). This term is 
analogous to the second term in square bracket in eq. (2.4o) and would coincide with it if w = —1/3 
and, namely, if 9 = —h/2. If this is true, apart of a different power distribution along the I axis and 



(2.2) 
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Figure 1: TE correlation spectra [Cj = 1(1 + 1)\CJ e \/Ait) for ACDM (the same in the three plots) and 
RP. Solid (dashed) lines refer to r = 0.14 (0.20). In the third plot, low-1 peaks are negative. 



geometric effects at greater /, there can be similarities in the behavior of open and dynamical DE 
models. The relation between h/2 and Ode can then be studied through the equation 



Ode + \ = -— I; — [Sde + 3-(cg - w)5 DE ], 



(2.5) 



whose validity indicates that DE behaves as a fluid. An order of magnitude estimate, however, tells 
us soon that Ode ~ k 2 t5DE', then, the ratio between Ode and the r.h.s. is ~ (t/L) 2 , where L is the 
scale related to k. Before horizon crossing (t <C L), the Ode term is negligible, in comparison with 
the r.h.s.. Hence, h/2 equates the r.h.s. and, therefore, the ratio —Ii/20de exceeds unity. At horizon 
crossing such ratio must approach unity and keep such as t grows greater than L. The main differences 
between open and dynamical DE models, in the r.h.s. of eqs. (2.4) are therefore relegated to times 
before horizon crossing. Afterward, the residual difference is due to a factor 1 + w. 

The changes in f/ directly act on F 2 . The equation fulfilled by this spectral component reads: 



F 2 (k,t) 



-F 2 {k,t) 



U(k,t) 
10 



+ |[2Fi(M) 



3F 3 (M)] 



fj(k,t) 



h(k,t) 



(2-6/) 



in flat models, while in open models we have: 

U(J3,t) 



F 2 (J3,t) = -f 



-F 2 (P,t) + 



10 



+ £[2b 2 F 1 (0,t)-M 3 F 3 (P,t)] + ^b 2 U(/3,t) + ^Mj . (2 . 6o ) 



Accordingly, when open models show negative TE correlations at low I, we expect something similar 
in flat RP models. On the contrary, models like SUGRA, with a cosmic acceleration closer to ACDM 
models, are not expected to give negative TE correlation. 

Let us then discuss how the likelihood distribution can be obtained from model spectra, taking 
into account that the number of pixels for anisotropy and polarization (Nt and Np) in our (artificial) 
data are different. Let Tj be the anisotropy data measured in Nt pixels and Qj and Uj be the Stokes 

parameters in Np pixels. In general, let be x = (T\ , , T/v T , Q\ , ,Qn p ,Ui, , Un p ) • x is a vector 

of N s = Nt + 2Np components, defining an observed state of anisotropy and polarization. Once a 
model is assigned, the Cj are uniquely determined. Passing to a data vector x, instead, amounts to 
performing a model realization. Vice-versa, once the N s component data vector d is given, the model 
is not uniquely fixed. 
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Figure 2: Histograms give the distribution of peak likelihood. Continuous curves give the likelihood 
distribution averaged over realizations (arbitrary but equal normalization). Pixel noise is 0.2 //K in 
the upper panels and 2 /j,K in the lower panels; all plots refer to r = 0.14. 

Then, the likelihood of a model, whose angular spectra are Ci, when the data d are observed, reads 

L{d\Cf ) oc [det M]~5 exp [ — ^d T M x d] . (2.7) 

The main ingredient of L is the correlation matrix M;j = (xf Xj) = Sy + N^; here Sjj is the signal 
term and Njj is due to the noise. The components Mij yield the expected correlation between the ith 
and jth elements of data vectors x corresponding to particular choices of C\ [10]. The construction 
of the noise term is simpler, as we expect no noise correlation, and the matrix Njj = Sij (JT,pix (for 
i = 1, Nt) and = 5^ ap, p i x 2 (for i = Nt + 1, N s ) is diagonal. 

3 Angular spectra and likelihood: results. 

In Fig. 1 we show the spectra Cf E for ACDM and RP models, if r = 0.14 and r = 0.20 (WMAPfll] 
suggests that r = 0.17 ± 0.04). The value of A = log 10 (A/GeV) is given in top of the frames. Low-/ 
anticorrelation is found only for A >~ 10 and is shown in the third plot. The differences between 
ACDM and RP increase for growing A, but are already significant even for A = 2. 

At low I, cosmic variance must be taken into account, aside of instrumental variance. We then 
perform a large number (1000) of realizations of RP sky models for A =2, 5 and 8, for r = 0.14 and 
0.20. In Fig. 2 (3) we report report results for r = 0.14 (0.20). The values of A are shown in top of 
each plot. Lower and upper plots correspond to a = 2 and 0.2 /xK, respectively. The former value 
approaches WMAP and SPOrt expectations; the latter value might be approached by the PLANCK 
[12] experiment. 



4 Conclusions 

These figures show that, as expected, RP models are more easily distinguishable from ACDM for lower 
a and higher A. We can assume that a RP model gives a signal different from ACDM when A > is 
detected. However, even in the less favorable case considered, when a = 2/xK, r = 0.14 and A = 2, 
the peak likelihood is at A > in ~ 72 % of cases. For A = 5, this fraction reaches ~ 91 %. Likelihood 
distributions, averaged both over cosmic and instrimental variances, tell us that RP models begin to 
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Figure 3: The same as Fig. 2, but for an optical depth of r = 0.20. 

be distinguishable from ACDM, at 1-a level, for A ~ 5 and r = 0.20, if a = 2fiK. For a = .2/j,K, 
instead, this is already true for A = 2. See the figures for further details. 

A general conclusion that can be drawn from this analysis is that cosmic variance is far from 
being a serious limit to model parameter detection from large angle spectral analysis, at the present 
sensitivity levels. Even for a sensitivity improved by a factor 10, hystograms are still contained inside 
the curve, showing that there is a relevant space for further improvements. 
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